Liquid drop splashing on smooth, rough and textured surfaces 
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Splashing occurs when a liquid drop hits a dry solid surface at high velocity. This paper reports 
experimental studies of how the splash depends on the roughness and the texture of the surfaces 
as well as the viscosity of the liquid. For smooth surfaces, there is a "corona" splash caused by 
the presence of air surrounding the drop. There are several regimes that occur as the velocity and 
liquid viscosity are varied. There is also a "prompt" splash that depends on the roughness and 
texture of the surfaces. A measurement of the size distribution of the ejected droplets is sensitive 
to the surface roughness. For a textured surface in which pillars are arranged in a square lattice, 
experiment shows that the splashing has a four-fold symmetry. The splash occurs predominantly 
along the diagonal directions. In this geometry, two factors affect splashing the most: the pillar 
height and spacing between pillars. 

PACS numbers: 47.20.Cq, 47.20.Ma, 47.40.Nm, 47.55.D- 
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I. Introduction 

When a liquid drop hits a solid surface, it often 
splashes and breaks into thousands of smaller droplets. 
Splashing is an excellent example of a singular breakup 
phenomenon with an underlying instability that is still 
not properly understood. As illustration, it was only 
recently discovered that the surrounding air pressure 
is an important parameter for creating a splash on 
a smooth dry substrate so that the splash can be 
completely suppressed in a low pressure environment 
1] . Splashing is also broadly important in industry with 
applications in ink-jet printing[2|, combustion of liquid 
fuel[||], spray drying^ and surface coatingQ. 

There are two distinct types of splashing "corona" 
and "prompt" . Corona splashing occurs on smooth sur- 
faces, where a symmetric corona is first formed, and 
droplets are ejected from the expanding corona; prompt 
splashing takes place on rough surfaces, where there is no 
corona, and droplets are created at the spreading contact 
line. Fig. 1 shows photographs of the two cases. A pre- 
vious study proposed to explain this difference: corona 
splashing is caused by the effects of the air surrounding 
the drop and prompt splashing is caused by the effects of 
surface roughness 

Splashing has been studied since the time of Worthing- 
ton in 1867 [8[. Since then, there have been many ex- 
perimental studies which have attempted to find a cri- 
terion for when splashing would occur. Notably, Mundo 
et a?[9( established an empirical relationship for the no- 
splashing to splashing transition that depended on the 
surface roughness, R a , the velocity of impact, Vq, the 
surface tension of the fluid, a, the diameter of the drop, 
D, the dynamic viscosity of the fluid, /i, and the density 
of the fluid, p. They found that: We 1 l 2 -Re 1 ^ = K c [R a ], 
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FIG. 1: Corona splash and prompt splash. The photo- 
graph on the left is a corona splash on a smooth dry sur- 
face. Droplets are created from a symmetric corona. In the 
photograph on the right, a prompt splash occurs on a rough 
dry surface. In this case there is no corona and droplets are 
ejected from the advancing contact line. 



where We and Re are Weber number and Reynolds num- 
ber respectively: We = pDV^/a, Re = pDVo/p. K c is 
a constant and that depends on surface roughness R a . 
Splashing will occur when K > K c ~ 50. When K < K c 
no splashing will occur. Wuflfi| and Range et a/[ll[ in- 
vestigated the dependence of splashing on the Ohnesorge 
number, Oh = pj \JDop. They studied the Oh << 1 
case , where they could neglect the effects of viscosity, 
and obtained the relationship: We c — a ■ log b R a , where 
a and b are fitting parameters. When We > We c , they 
saw a splash. Note that neither of these relationships 
take into account the effects of the gas surrounding the 
liquid during the splash. 

Some researchers also investigated the fingering insta- 
bility at the rim of the expanding liquid disc. Allen [l2| 
proposed that the Rayleigh- Taylor instability caused the 
fingering. Bhola et al [l3j ] and Mehdizadeh et al [3] ob- 
tained reasonable agreement between this theory and 
their experiments. Thoroddsen et al\lE^ experimentally 
studied the fingers and proposed that the instability is 
caused by the presence of air trapped under the liquid 
drop. 

This paper reports on experiments both for corona 
splashing on smooth dry surfaces and for prompt splash- 
ing on rough and textured dry substrates. For corona 
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splashing, there are several regimes that depend on the 
velocity of impact and the fluid viscosity. Undulations 
around the rim of the spreading fluid are measured as a 
function of the air pressure on the smooth dry surfaces. 
There is a sharp jump in the number of undulations 
at the threshold pressure. For prompt splashing, both 
random roughness and roughness created by regular 
textured surface were studied. For a textured surface 
consisting of a regular array of pillars, the dependence 
of splashing on the vertical pillar height, lateral pillar 
size and pillar spacing was studied independently. 

II. Corona splash on smooth surface 

Previous experiment has shown that the sur- 
rounding air is crucial for corona splashing on a smooth 
dry surface [1]: The rows of Fig. 2 show images of 
splash at different background air pressures for a drop 
of ethanol hitting a glass substrate. Surprisingly, as the 



pressure is lowered, fewer droplets are ejected; under low 
enough pressure no droplets emerge at all after impact. 
At a threshold pressure, Pt, the splash just begins to be 
formed as is shown in the second row of the figure. 

The threshold pressure, Pt, as a function of impact 
velocity, Vo, is shown in the main panel of Fig. 3. The 
curve is not monotonic. In the high velocity region above 
a characteristic velocity, V* , Pt decreases as the impact 
velocity is raised. This is what we might naively ex- 
pect. However, in the region Vq < V*, the curve is non- 
monotonic. This non-monotonicity indicates two differ- 
ent regimes at low and high velocities. Further experi- 
ments show that V* varies with liquid viscosity and drop 
size [16} . 

Experiments have also revealed that when the sur- 
rounding gas is heavier (for example, using Kr and SFg) 
and with larger liquid viscosity it is easier to create a 
splash. We compared two stresses the destabilizing 
stress from air, Eg, and the stabilizing stress from sur- 
face tension, Ex,, and found: 
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FIG. 2: Photographs of a liquid drop hitting a smooth dry 
substrate. A 3.4±0.1 mm diameter alcohol drop hits a smooth 
glass substrate with impact velocity Vo = 3.74 ± 0.02 m/s at 
different background air pressures. Each row shows the drop 
at two times, 0.276 ms and 0.552 ms after impact. In the top 
row, P = 100 kPa (atmospheric pressure), the drop splashes. 
In the second row, at the threshold pressure, Pt = 38.4 kPa, 
the drop emits only a few droplets, traveling at a small angle 
with respect to the surface. In the third row, at P = 30.0 
kPa, there is no splashing but there are undulations at the 
rim. In the fourth row, at P = 17.2 kPa, there is no splashing 
and no apparent undulations in the rim of the drop. Taken 
from 0]. 
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Here 7 is the adiabatic constant of the gas, Mq is the 
gas molecular weight, ks is Boltzmann's constant, T is 
the temperature, D is the diameter of the drop, vt is 
the kinematic viscosity of liquid, and a is the surface 
tension. A heavier gas or a larger liquid viscosity will 
increase the ratio, Eg/El . The ratio of these two stresses 
was found to be approximately constant for velocities 



CO 
CL 




V o (m/s) 

FIG. 3: Threshold pressure versus impact velocity and the 
collapse of data. Main panel shows Pt vs. Vq in air. The 
curve is nonmonotonic: seemingly two distinct regimes are 
separated by a velocity, V* . The inset plots Eg/£l versus 
Vo at threshold pressure Pt, in the region Vo > V* , for gases 
of different molecular weight, Mg'- He(Mir e = 4), air(M a ; r = 
29), Kx(M K r = 84), SF 6 (M s f 6 = 146) and for liquids of differ- 
ent viscosity, vl: Methanol(i^Meih = 0.68cSt), Ethanol(i/E to h 
= 1.36cSt), 2-Propanol(^2-Pro = 2.60cSt). At threshold pres- 
sure, all Eg/S_l collapse approximately onto a constant value, 
0.45. Taken from Q. 
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above V* at threshold pressure. This is shown in the 
inset to Fig. 3 where Yiq/Hl at threshold pressure Pt 
is plotted for gases of different molecular weights (4 ~ 
146 Dalton), liquids of different viscosities (0.68 ~ 2.6 
cSt) and different impact velocities (2.5 ~ 7 m/s). At 
threshold pressure, in the regime V > V* , Sg/Sl = 0.45 
so that Eq.l successfully collapses all the data without 
any fitting parameter. 

The prediction of Eq.l that increasing the liquid vis- 
cosity leads to a lowering of the threshold pressure was 
verified by the data in Fig. 3 which spanned the range 
0.68cSt <v L < 2.60cSt. Nevertheless this result is coun- 
terintuitive from our experience with high viscosity liq- 
uids and calls for more experiments covering a broader 
range of viscosity. By using silicone oils of different 
molecular weights, the liquid viscosity could be varied 




by more than one order of magnitude, while keeping a 
very similar mass density (0.82 ~ 0.95g/cm 3 ) and surface 
tension(17.4 ~ 21mN/m). Fig. 4 shows photographs of 
a relatively viscous (5 cSt) silicone oil drop hitting a dry 
glass substrate under different pressures of air. Again, we 
find that the splash decreases as the air pressure is de- 
creased, and that no splash occurs when the background 
pressure is low enough. But one difference between Fig. 4 
and Fig. 2 is that splashing occurs at a much later time 
when the viscosity is large. This is most obvious if one 
compares the second rows. Clearly one effect of viscosity 
is to delay the splashing time, as we might have expected. 

Fig. 5 shows the threshold pressure, Pt, vs. liquid vis- 
cosity, vl, for 3.1 ± 0.1 mm diameter drops hitting the 
substrate with an impact velocity Vq — 4.03 ± 0.05 m/s. 
The upper curve shows, as before, the splashing thresh- 
old pressure, Pt, where splashing is first detected. The 
lower curve shows the threshold pressure, Pr—bump, for 
where an undulation in the expanding sheet of liquid is 
first observed. Pr-bump is defined as the lowest pressure 
at which undulations (or bumps) first show up, and be- 
low which no undulations can be seen. Both threshold 
pressures first decrease then increase with increasing vis- 
cosity. This indicates two different regimes. At low v L , 
as the viscosity is increased, the threshold pressure to 
create a splash decreases. Thus viscosity helps to pro- 
duce splash, as predicted by Eq.l. The solid line is the 
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FIG. 4: Splash of a viscous drop(z/ L = 5cSt). A 3.1 ± 0.1 
mm diameter silicone oil drop hits a smooth glass substrate 
at impact velocity Vb = 4.03 ± 0.05 m/s under different air 
pressures. Each row shows the drop at two times: 0.735 ms 
and 2.478 ms after impact. In the top row, at P = 100 kPa, 
there is a pronounced splash. In the second row, at threshold 
pressure, Pt = 36.0 kPa, the drop just starts to splash. In 
the third row, at P — 30.7 kPa, there is no splash but there 
are undulations in the thickness of the rim. In the fourth 
row, P = 26.7 kPa, there is no splashing and no apparent 
undulations in the rim. The general property that less air 
leads to less splashing is similar to the low viscosity case shown 
in Fig. 2. However, at high viscosity splashing occurs at a 
later time (2.478 ms for the second row) than it does for low 
viscosity (0.552 ms for Fig. 2 the second row). 
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FIG. 5: Threshold pressures versus viscosity. Except for the 
point at the smallest viscosity, which is for methanol, the flu- 
ids were silicone oils with different molecular weights. Two 
curves are shown: (i) the splash threshold pressure, Pt(»), 
for a splash to appear and (ii) the bump threshold pressure, 
PT-bump(°), where the bump at the leading edge first ap- 
pears. The curves are both non-monotonic. At low viscosity, 
the threshold pressures decrease with increasing ul whereas 
at high viscosity, they increase with ul- The solid line is the 
curve predicted by Eq.l: Pt ~ 1/^/vl. The curve fits the 
small viscosity regime very well, but does not capture at all 
the trend at high vl- In these experiments, the impact veloc- 
ity and drop diameter are kept fixed at Vb = 4.03 ±0.05 m/s, 
D = 3.1 ±0.1 mm. 
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scaling relation derived from Eq. 1 , which agrees well with 
the small i>i, data. However, the prediction starts to de- 
viate at higher viscosity where the threshold pressures 
increase with ul- In this regime, the higher the viscosity, 
the higher the pressure of air needed to create a splash 
and viscosity suppresses splashing. 

Why there are two different behaviors? We think 
that for the low vl regime, the expanding liquid film 
is stabilized mainly by surface tension so that viscosity 
only affects the film thickness: d ~ \fvj^. Thus a larger 
Vjj causes a thicker film which is easier to destabilize. 
But for the high vl regime, viscous drag is important 
and helps to stabilize the spreading drop. 

III. Number of undulations versus pressure 

Researchers have extensively studied the finger- 
ing instability that occurs as a splash is created 
EH [3 EH • "Fingers" mean long protrusions at the 
rim of the expanding liquid film. In our experiment, we 
also observed undulations around the rim as shown in 
the inset to Fig. 6. Here, a feature similar to "fingering" 
is observed. However in this case, the undulations 
do not extend very far out from the rim. In order to 
prevent possible confusion, we call them "undulations" 
or "bumps" . Previous studies have focused on the 
number of fingers as a function of impact velocity and 
surface roughness. Here we concentrate on the behavior 
and number of the bumps as the air pressure is varied. 
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FIG. 6: Number of bumps as a function of air pressure. Inset 
is a bottom-view photograph showing the undulations. The 
number of undulations shown in the main panel, Nt ump , is 
counted from such images. There is a sudden change in Nbump 
around Pr-bump = 29 kPa. At that pressure Nbump jumps 
from to a finite value and stays constant above that pressure. 



Nbump jumps from zero to a finite value, and stays con- 
stant for higher pressures. The absence of undulations 
at low pressure suggests that no instability can grow 
below Px-bump- The apparent pressure independence 
of Nbump above Pr-bump could be due to the narrow 
pressure range we are able to measure: above a certain 
pressure the entire expanding film is lifted into air. 

IV. Discussion on instability mechanism 

What is the mechanism for destabilizing the system 
and causing the occurrence of a splash on a smooth sur- 
face? What is the instability that eventually produces a 
splash? This question is still under debate. One prevail- 
ing theory first proposed by Allen 12], is that it is due to 
the Rayleigh- Taylor instabilitv[17|. This is an instabil- 
ity occurring at the interface of two fluids with different 
densities, when the acceleration, a, points from the light 
fluid, p\, towards the heavy fluid, p 2 . Linear instabil- 
ity theory predicts an exponential growth in amplitude, 
A ~ Exp(ct), with growth rate c. The wavenumber, k m , 
and the growth rate, c m , of the fastest growing mode are 
predicted to be[Hj|: 



h — 



a(p 2 - pi) 
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I 2ak m (p 2 - pi) 
3(P2 + Pi) 



(2) 
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Here the light fluid is the surrounding air and the heavy 
fluid is the expanding liquid. The expanding disc de- 
celerates so that the direction of a satisfies the require- 
ment for the Rayleigh- Taylor instability. However, since 
Pi » Pi, Eqs.2 and 3 do not vary appreciably with air 
pressure. Therefore splashing should not vary with air 
pressure if it were caused by Rayleigh- Taylor instability. 
This is inconsistent with our experiment. 

Another interface instability, the Kelvin-Hclmholtz in- 
stability, can take place when there is velocity jump at 
the interface. For inviscid fluids and p 2 » pi, the 
wavenumber and the growth rate of the fastest growing 
mode are: Q1, [11 SIl : 



_ 2 Pl u 2 
3 a 



(4) 



In Fig. 6 we show Nbump at different pressures. We 
determine the number of undulations, Nbump, from 
pictures such as the one shown in the inset. At all 
pressures, Nbump is measured at the same fixed radius 
of expansion where the undulations are most clear. The 
main panel shows that at threshold pressure Pr-bump, 



(5) 



With u the relative veloci ty betwe en two fluids at the 
interface. In our case, u ~ y DV$/At is the velocity of the 
expanding liquid film. The Kelvin-Helmholtz instability 
strongly depends on the density of the lighter fluid, pi, 
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and thus may be relevant to our experiment. However, 
our previous results [l| indicate that compressibility of air 
is important. This suggests that we should replace the 
Bernoulli term p±u 2 in Eq.4 with piCcu, with Co — 
y/jkBT/Mc the speed of sound in the surrounding gas. 

2 Pl C G u 2P RM^ [DVo 

fcm i— = 37 yi^r y— (6) 

The characteristic length in the expanding liquid film 
is the film thickness, d. This suggests that the instability 
might be able to grow if: 

k m ~ 1/d (7) 

Since we also have d ~ yJviX, from Eq.6 and Eq.7, we 
obtain as a criterion for the instability to grow: 



this behavior in the following manner. After impact, 
the thickness of expanding film, d, grows continuously 
from being molecularly thin just after the impact, 
to approximately 40 ~ 50/xm at the end of the film 
expansion. When the surface roughness is small, d can 
grow to be much larger than R a . At the beginning, when 
d is small, the film is thinner than R a and continues 
to eject droplets until the film becomes much thicker 
than the surface roughness. After d grows to be larger 
than R a , the roughness is too small to destabilize the 
liquid film and produce a splash. Thus the distribution 
of ejected droplets reflects the surface roughness R a and 
we find ro rs R a - However, when the roughness is large, 
d can never grow to be greater than R a . Consequently, 
ro can only grow to the maximum size of d at its final 
thickness. This is consistent with the decay constant, 
ro, saturating around 40/zm, which is roughly the film 
thickness at the end of expansion. 



DVo 
Ak B T 



1 



(8) 



The left hand side is exactly Eg/El in Eq.l. In the 
low viscosity regime, our experiment gives Eg/El = 0.45 
for the splashing threshold which is consistent with the 
criterion in Eq.8. This suggests the possibility the 
Kelvin-Helmholtz instability may be the underlying 
instability mechanism for corona splashing. 

V. Prompt splash on rough surface 

A completely different type of splash, the prompt 
splash, occurs on rough surfaces. By systematically 
varying the degree of surface roughness and the air 
pressure, we discovered two different mechanisms for 
the two kinds of splashes: surrounding air is responsible 
for the corona splash discussed above and surface 
roughness is responsible for the prompt splash Q. Under 
ordinary conditions (atmospheric pressure and non-zero 
roughness), a splash is a mixture of both contributions. 
By working under low pressure with a negligible amount 
of air, we are able to study pure prompt splashing. 

Since a prompt splash is caused by surface roughness, 
it may retain information about surface roughness in the 
distribution of sizes of the ejected droplets. We find 
that this is the case. We mix a small amount of ink 
into our ethanol and then collect the ejected droplets 
on a sheet of white paper. We then obtained the sizes 
of ejected droplets by measuring the size and darkness 
of the stains left on the paper. Our previous study[7J 
shows that the number of droplets, N, decays exponen- 
tially with their radius, r: N(r) ~ Exp{— r/ro)(see inset 
of Fig. 9(b)). This indicates the existence of a charac- 
teristic decay length, tq. 

The decaylength, ro, correlates with roughness of the 
surface, R a '7]. At small R a , we have the relationship 

~ Ra', but for large roughness, this breaks down as 
ro saturates at a constant value. We can understand 



VI. Prompt splash on textured surface 

The last section showed that roughness has a 
strong effect on prompt splash. To understand this 
dependence in more detail and to understand how 
surface properties affect splashing, we study splashing 
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FIG. 7: Prompt splashing on a textured surface. (a)Top row 
shows a side view of prompt splashing on a textured surface. 
Bottom row shows a bottom view of prompt splashing on a 
textured surface. There is a clear four-fold symmetry in the 
splash which is predominantly in the diagonal directions of the 
square lattice created by the pillars, (b) Making a textured 
surface with UV-lithography technique. Left cartoon shows 
the UV-lithography process. The right picture shows a typical 
textured surface under the microscope. We define the pillar 
height as h, lateral pillar size as I and spacing between pillars 
as s. For this particular substrate, h = 18pm, I = s = 60/im. 
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on a well-defined textured surface of regular patterns. 

The textured surface is made with UV-lithography: 
We first spin coat UV epoxy (SU8-2000, MicroChem 
Corp.) onto a clean glass microscope slide. We then cover 
the slide with a mask of pre-designed pattern (square 
blocks in a square lattice) and expose the slide to UV 
light. After development, the UV epoxy film which is di- 
rectly under the transparent part of the mask will harden 
and the rest of it can be rinsed away, resulting in a struc- 
ture on the substrate of square pillars arranged in a 2D 
square lattice as shown in Fig. 7b. There are three im- 
portant quantities relating to our splash experiments in 
this textured surface: (f ) the vertical pillar height, h, (2) 
the lateral pillar size, I, and (3) the lateral spacing be- 
tween pillars, s. By changing the spin speed, we can vary 
pillar height, h; by designing different mask patterns, we 
can vary both I and s independently. Thus we can vary 
every aspect of the structure. 

Fig. 7(a) shows photographs of a prompt splash on 
a textured surface under low air pressure. The impact 
velocity is 4.3 ± 0.1m/ s, and the drop diameter is D = 
3.4±0.1ram. The top row shows a side view of the splash. 
It has a similar look as the prompt splash on an ordinary 
rough surface. However, the bottom view shown in the 
second row reveals a very striking feature: the splashing 
occurs with four fold symmetry. The droplets are ejected 
predominantly along the diagonal directions of the square 
lattice. Fig. 7(b) shows the process of UV-lithography 
and a picture of textured surface under a microscope. 

We can now vary the profiles of the surface and de- 
termine their effect on the splash and the ejected droplet 
distribution. Again we use the ink spot technique to mea- 
sure the size distribution of droplets, as mentioned in last 
section. In the first set of experiments we first keep lat- 
eral size constant at I = s = 60/^m, and vary the vertical 
height of the pillars, h. Fig. 8(a) shows the number of 
droplets, N, versus their radius, r. Similar to the case 
with random roughness, we find an exponential decay at 
large r, N(r) <~ exp(—r/r ), with a characteristic decay 
length, ro, that varies with h. Fig. 8(b) shows that ro 
varies with ft, in a nonmonotonic manner. For small h, ro 
increases with, and has a value comparable to, h. This 
indicates that ro is determined by h. However, when h 
is greater than 18/zm, the opposite trend occurs: ro de- 
creases as h increases. Fig. 8(c) shows the sum of the 
areas created by all of the ink spots, A tot , as a function 
of h. A to t is a quantity that indicates the total amount of 
ejected droplets. Fig. 8(b) and (c) have the same shape, 
indicating that r and A tot are strongly correlated. The 
decreasing trend for large h implies that larger roughness 
leads to less splashing. When h is greater than 60/im, 
there is no splash at all. 

These results are counterintuitive. We suspect they are 
caused by the way in which the impacting liquid drop can 
flow between the channels set up by the pillar structure. 
At small pillar heights, the liquid in the drop can eas- 
ily reach the bottom of the canyon between the pillars 
and can then expand along the bottom surface. During 
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FIG. 8: Decay length and total splash area versus pillar 
height, (a) The number of ejected droplets, N, as a function 
of droplet radius, r, for three pillar heights, h. The lateral di- 
mensions are kept fixed at I = s = 60/im. The straight lines 
are an exponential fit to the large r tail of the distribution: 
N ~ Exp(-r/r ), with r = 0.0017mm ( — ), 0.0098mm (— 
) and 0.005mm ( — ) for h= 6/ttm(A), h = 30/im(»), and h 
= 54/um(D) respectively. Note that ro does not vary mono- 
tonically with h. (b) The exponential decay length, ro, versus 
pillar height, h. ro first increases then decreases as h increases, 
(c) The total ink spot area, A to t, as a function of h, for the 
same set of experiments shown in (b). The curve in (c) has a 
similar shape as (b) indicating that A to t and ro are strongly 
correlated. 



expansion, the liquid film is destabilized by the pillars, 
producing droplets with a size related to the pillar height, 
h. This produces a positive correlation between ro and 
h for h < 18^m. However, as h increases, it is increas- 
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ingly difficult for the impacting drop to reach the bottom 
surface. In this case, we suspect that much of the drop 
expands on top of the pillars, rather than between them. 
Once h sw 60/im, that is when the height is about the 
same size as the lateral dimensions I and s, the situation 
resembles a drop expanding on a flat surface with many 
holes rather than pillars. Here all the obstructions are 
underneath the liquid film and make only a small pertur- 
bation to its expansion. Because the impacting drop can 
only penetrate a finite depth below the pillar top, it does 
not know how far away it is from the bottom surface. 
This suggests that the amount of splashing should satu- 
rate as the pillar height is increased. We do not have a 
good explanation of the surprising fact that the splashing 
can be completely eliminated if the pillars are sufficiently 
tall. 

To some extent, this is similar to the Cassie state of 
a drop on a superhydrophobic rough surface studied by 
Quere et al [23 . |23| , where a water drop can sit on top of 
air trapped in the rough profile of the substrate. However 
we note that their case is static whereas ours is proba- 
bly driven by the fast dynamics of the expanding drop. 
Moreover in the case studied by Quere et al the air plays 
an important role in supporting the weight of the drop. 
In our situation air has been pumped out of the system. 

To understand the effect of the lateral dimension on the 
splashing, we make substrates of different / and s, while 
keeping h fixed. Fig. 9(a) main panel shows iV versus r 
with an exponential fitting function, N ~ Exp(—r/ro), 
for different lateral sizes. Fig. 9(b) plots ro as a function 
of lateral pillar size I and spacing s. As I and s are varied 
we keep I = s, and h = 10^m. Fig. 9(c) shows the total 
area, Atot, vs. I and s. A to t has the same dependence 
on I and s as does r$. Both quantities increase with 
lateral size in most of our range, then decrease at the end. 
This means that increasing the lateral dimensions will 
enhance splashing for small I and s. When the pillars are 
too sparse, splashing becomes less pronounced suggesting 
that it is more difficult to destabilize the liquid film. We 
should also note that ro is always much smaller than I 
and s, while much closer to the pillar height h = lO/im. 
This indicates that h is more important in determining 
ro than are I and s. 

A comparison of a textured surface with a random 
roughness surface is shown in the inset of Fig. 9 (a). Both 
curves decay exponentially, but random roughness curve 
(upper curve) has a much larger ro. This is surpris- 
ing because both curves have similar roughness (I — 
s = 20/ito, h — 10/im for the textured surface and 
R a = 16/ito for the case of random roughness). We can 
understand this qualitatively using the data shown in 



Fig. 8(b). There we see ro starts to decrease at h 
Thus here at h = 10 (im = \l, it is already difficult for 
the drop to reach the bottom of the substrate and get 
destabilized. However, the random roughness is made 
by particles coated on surface. This substrate never re- 
sembles a flat surface with many holes. Therefore the 
random roughness surface makes a much larger splash 
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FIG. 9: Decay length and total splash area versus lateral 
size of the pillars, (a) TV as a function of r for three lat- 
eral sizes. The pillar height h is fixed at h — W^im. The 
straight lines are an exponential fit to the large r tail of the 
distribution: TV ~ Exp(—r/ro), with ro = 0.0014mm ( — ), 
0.010mm ( — ) and 0.0061mm ( — ) for I = s = 20/im(A), 
1 = 8 = 80^m(»), and I = s = lOO^m(D) respectively. 
Note that ro does not vary monotonically with I and s. In- 
set compares the distribution for a textured surface with 
I = s = 20/im, h — 10/im(A) with a sample with comparable 
random roughness R a = 16/im(o). The random roughness 
creates a much larger splash, (b) ro is plotted versus I and s. 
ro first increases then decreases with / and s. (c) Atot versus 
I and s. Again we see a similar shape as in (b). 



with a larger ro. 

Fig. 9 demonstrates that ro changes as we vary I and s 
together. One further question is whether this is caused 
by a change in / or a change in s or in both? We can check 
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FIG. 10: Effect on splash characteristics of varying I as com- 
pared to the effect of varying s. (a) We independently vary 
lateral pillar size , l(x), and spacing between pillars, s(o), to 
compare their effect. Two lines are guides to the eye for I ( — 
-) and s ( — ). Pillar height is kept fixed at ft = 10/im. When 
I is varied, s — 60fim is held fixed; when s is varied, I = 60/im 
is held fixed. The comparison shows that ro changes with s 
but not I. (b) A plot of the total area, A to t, as a function of I 
or s, for the same experiment as in (a). It produces the same 
trend as in (a). 



this by varying I and s independently. Fig. 10 shows the 
result. The different symbols are for varying I and s 
separately while leaving all other conditions unchanged. 
Apparently increasing the spacing between pillars, s, en- 
hances ro and A to t while I has a much smaller effect on 
the splashing behavior. 

The fact that increasing s enhances splashing helps to 
explain why we see splashing in the diagonal directions in 
Fig. 7a. Along the diagonal, the distance between pillars 
is the greatest. Because s is largest in those directions, 
splashing preferentially occurs in those directions. 

The textured substrate not only affects the prompt 
splashing caused by surface roughness, it also changes 
the behavior of corona splashing caused by the surround- 
ing air. Fig. 11(a) shows a typical corona splash on a 
smooth surface at atmospheric pressure, while Fig. 11(b) 
shows, at the same pressure, no splash at all on a textured 
surface consisting of tall pillars. In both experiments 
the drop hits the substrate at the same impact velocity. 
Different amounts of splashing can also be achieved by 
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FIG. 11: Effect of pillars on suppressing corona splashing. 
All experiments are done at atmospheric pressure, with Vo = 
4.3m/s and D — 3Amm. (a) A corona splash on a smooth 
surface, (b) Splashing is completely suppressed on a surface 
consisting of pillars with I = s = 60pm and h = 125^im. We 
can tune the amount of splashing by varying the pillar height 
even under atmospheric pressure. 



creating pillars with intermediate heights. These results 
suggest that the pillars form channels through which the 
air can escape so that the importance of the air for creat- 
ing the splash is minimized. This discovery demonstrates 
another way in which one can suppress splashing. More- 
over, it has the advantage that this suppression can be 
achieved without decreasing the gas pressure. 

VII. Conclusion 

This paper systematically studied the splashing of liq- 
uid drops on various dry solid surfaces. This study cor- 
roborates that there are two mechanisms corresponding 
to the two kinds of splashes. Air causes the corona splash 
on smooth dry surfaces and substrate roughness causes 
the prompt splash. For the corona splash, we discovered 
several regimes. At high impact velocity, there are two 
regimes as the viscosity of the liquid is varied. We also 
studied the finger instability as function of air pressure 
and find a jump in the number of bumps. We suspect that 
Kelvin-Hclmholtz instability coupled with the compress- 
ibility of air, is a possible mechanism for the splashing 
instability. This mechanism agrees well with our experi- 
mental data. 

In order to examine the effect of surface roughness, 
we studied splashing on textured surfaces consisting of 
square pillars arranged in a square lattice. We found 
that the dimensions of the pillars strongly affect splash- 
ing. Here the pillar height, h, is found to be the most 
important factor determining characteristic decay length, 
ro. We discovered that the splash preserves the symme- 
try of substrate. This shows that the splash direction 
can be controlled. We also find that corona splash un- 
der atmospheric pressure can be suppressed by making 
tall pillars on surface. This provide another way to re- 
duce splash even under normal pressure. Since splashing 
is involved in many industrial processes 0, [H, 0, Q, these 
discoveries could have important practical applications. 
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